Wigner Crystals in the lowest Landau level at low filling factors 
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We report on results of finite-size numerical studies of partially filled lowest Landau level, at low 
electron filling factors. We find convincing evidence suggesting that electrons form Wigner Crystals 
at sufficiently low filling factors, and the critical filling factor is ~ 1/7. At = 1/7 we find the 
system undergoes a phase transition from the Wigner Crystal to the incompressible Laughlin state 
when the short-range part of the Coulomb interaction is modified slightly. This transition is either 
continuous or very weakly first order. 
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Interacting electrons at sufficiently J ow densities and 
temperatures were shown by WignerEI to form a crys- 
talline state. This ground state results from the domi- 
nance of the potential over the kinetic energy which oc- 
curs in the dilute limit. A Wigner Crystal (WC) state 
has been observed in a low-density two- dimensional (2D) 
electroii gas system trapped on the surface of a helium 
liquidfl While unlikely to be observed at metallic densi- 
ties in 3D this crystalline phase is in fact expected to oc- 
cur in the 2D electron gas confined in semiconductor het- 
erostructures or quantum wells, and subject to a strong 
perpendicular magnetic field (such that all electrons are 
in the lowest Landau level (LLL)). Under these condi- 
tions the kinetic energy is completely quenched making 
it an ideal environment in which to observe Wigner crys- 
tallization; the WC state is expected to have the lowest 
Coulomb interaction energy when the Landau level filling 
factor {v) is sufficiently low. Indeed, substantial exper- 
imental evidence, mostly from transport measurements, 
points-.to a pinned WC state in the insulating phases of 
low vU However the situation is not yet definitive as di- 
rect experimental probe of the lattice structure of the 
WC is still lacking. In addition, the ubiquitous disorder 
in these systems could lead to other insulating phasesti 
such as the Hall insulator, or just complicate the detec- 
tion jsf the crystal itself. On the theoretical side, early 
workcl suggested that a 2D electron gas is always unsta- 
ble against formation of a WC when subject to a strong 
magnetic field, at zero temperature. However, the dis- 
covery of fractionaJ quantum Hall effect (FQHE) has led 
to the realizationd that electrons can also form incom- 
pressible liquid states, especially at the primary sequence 
V = 1/m, where "Vi= 3,5,---. By comparing the en- 
ergy of the Laughlind state with that of a correlated WC 
state (whose energy is considecably lower than that of a 
simple Hartree-Fock estimatecnZl) , Lam and GirvinH con- 
cluded that the critical filling factor I'c below which the 
WC forms is slightly above 1/7. This is in good agree- 
ment with transpuart experiments on clean samples. A 
more recent studya using a formalism based on composite 
fermionll^ (CF) wavefunctions calculated the collective 
excitation energy variationally and found that the ground 



state become unstable against proliferation of magneto- 
rotons (which points to formation of WC) for v < 1/9, 
although such an instability was not found at v = 1/7. 
A similar trend_.was observed earlier in the single mode 
approximation!^ where the magneto-roton gap while not 
zero was found to be diminished considerably at 1/7 and 
further at 1/9. This was interpreted as a precursor to 
magneto-roton mode softening and the transition to the 
WC phase. 

In this paper we present finite size exact diagonaliza- 
tion studies of systems with low v which directly reveal 
the formation of the WC for the first time in the LLL. We 
study systems with torus geometry. Our previous work 
on charge draisity wave (CDW) ground states in high 
Landau levela!30 has demonstrated that the torus geom- 
etry is advantageous as compared to, say, the spherical 
geometryliil, in studies of states with broken translational 
symmetry. The reasons are: (i) on a torus one can adjust 
the geometry to better accommodate a lattice in a finite- 
size system, (ii) The presence of translationaLsymmetry 
allows one to define a many-body momentumEj that can 
be used to label the eigenstates. By inspecting the mo- 
menta of the low-lying states, as well as the wavevector 
dependence of the density-density correlation functions, 
one can easily detect whether the translational symmetry 
is broken in the system, and determine lattice structure 
of the broken symmetry (crystal) state. 

We find strong evidence that WC forms at filling fac- 
tors < 1/7, and that the critical filling factor i'c is very 
close to 1/7. Furthermore, we have studied the com- 
petition between the WC state and the Laughlin state 
at z/ = 1/7 by varying the geometries of our finite-size 
systems, as well as the short-range part of the Coulomb 
interaction. We find that the Laughlin state is always 
continuously deformed into one of the low-lying states 
that correspond to the WC phase, when the geometry or 
the interaction is varied. This suggests that the quan- 
tum phase transition between the incompressible FQHE 
and compressible WC phases is either continuous or very 
weakly first-order. The methods used here are similar to 
our previous studies and provide a reliable way of detect- 
ing broken symmetry in finite systems. We first investi- 
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jate dependence of the spectrum on PBC geometry. 
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electrons, at filling factors v = 1/6,1/7, and 1/8, with 
rectangular geometry, for pure Coulomb interaction. For 
ly — 1/6 and 1/7 we plot energy levels for a series of 
aspect ratio oq, while for v — 1/8 we plot the energy level 
versus the magnitude of many-body momentum K, for 
Go — 0.35. One can see that for the = 1/8 case, as well 
a,s ly — 1/7 for a sizable range of ao: 0.3 < ao < 0.5, there 
are Nu = 6 nearly degenerate low-lying states (including 
the ground state) that form what we call the ground state 
manifold. The K's of the states in this manifold form a 
2D array, as we show in Fig. 2 for = 1/7 at oq = 3/8, 
indicating that the translational symmetry is broken in 
both directions, and the system has 2D crystalline order. 
These K's determine the lattice structure of the crystal. 
The number of unit cells included in the finite size system 

can be determined easilyEl: Nc = n'^ /Nd, where TV is 
the highest common divisor of N and the number of flux 
quanta N^. Here we have iV = 6 in both cases, thus 
Nc = 6 = TV, indicating there is one electron per unit 
cell. This is precisely what one expects for a Wigner 
Crystal (WC). This suggests that electrons form a WC 
at these filling factors, at least when a WC can be easily 
accommodated by the geometry of the finite size system. 
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FIG. 2. The allowed momenta K for the v = 1/7 system 
with N = & electrons with aspect ration ao = 3/8, and the 
K's of the ground state manifold. These low-lying K's form a 
nearly triangular array, which represent the reciprocal lattice 
of the Wigner Crystal that the electrons form. 



FIG. 1. The spectra of rectangular-shaped finite-size sys- 
tems with A'^ = 6 electrons, at filling factors 1/6, 1/7, and 1/8. 
The energy levels are plotted as functions of the aspect ratio 
for V = 1/Q and 1/7, and plotted as a function of the magni- 
tude of momentum K for v = 1/8, at aspect ratio 0.35. The 
degeneracy of each low-lying level is indicated by the integer 
beside it for z/ = 1/8. 

Fig. 1 shows the energy spectra of systems with iV = 6 



For V — 1/6, on the other hand, such a (near) degeneracy 
is absent for almost all geometries, except for perhaps at 
flo ~ 0.35. This suggests that the translational symmetry 
is not broken at v = 1/6, although the tendency toward 
WC formation is present, especially when the geometry 
of the finite-size system is favorable. We thus conclude 
that the critical filling factor Vc below which the Wigner 
Crystal forms is close to and probably slightly above 1/7, 
which is in excellent agreement with the estimate of Lam 
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and GirvinJa /^From the superlattice primitive basis vec- 
tors in the reciprocal spacetJ we obtain the parameters of 
the WC unit cell ai w 4.421, as ~ 6.655, and 9 « 131.6°. 
While this is not exactly a triangular lattice it is close 
to it. We believe the discrepancy is caused by the finite 
size of the system, and in the thermodynamic limit we 
expect to recover a triangular lattice. We emphasize the 
fact that the WC state is stabilized near oq = 0.4 is not 
an artifact of the geometry of the finite size system; as 
we will see below at filling factors = 1/3 and 1/5 where 
the system is an incompressible fluid and well described 
by the Laughlin state, no crystalline order is developed 
and the ground state has high overlap with the Laugh- 
lin state for all reasonable geometries, including the ones 
that are favorable for WC formation. 
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FIG. 3. The square of the overlap between the Laughlin 
state and the exact ground state of systems with rectangular 
geometry and a variety of filling factors and number of elec- 
trons, as a function of aspect ratio. The points at 1.2 are for 
hexagonal unit cells. 

At J/ = 1/m with m being an odd integer, it is knownEl 
that for certain short-range repulsive interactions, the 
Laughlin state is the exact ground state that describes 
an incompressible fractional quantum Hall liquid. In Fig. 
3 we show the square of the overlap between the Laugh- 
lin state with the exact ground state of the system for 
TO = 3, 5 and 7, at different aspect ratios and system 
sizes. It is seen that for = 1/3 and 1/5, the overlap 
is rather insensitive to the geometry of the unit cell and 
is close to 1, while for ly — 1/7 the square of the over- 
lap varies appreciably with the geometry, and dips below 
0.75 and 0.5 in certain range of gq for = 4 and = 6 
respectively. The sensitivity of the energy spectrum as 
well as the overlaps on the system geometry is a strong in- 
dication that the ground state is compressible. It reflects 
the fact that the system wants to form a WC, but is frus- 
trated when the aspect ratio oq is far from the optimal 
one for forming a crystal, ag ~ 0.4 for A^ = 6. It should 



be noted that the Laughlin state has K = 0, which is one 
of the K's in the ground state manifold, and in the case of 
a pure Coulomb interaction the ground state turn out to 
be in the K = sector for all the aspect ratios and sizes 
that we explored; thus the square of the overlap evolves 
continuously from being close to 1 near ao — 1, where the 
WC state is frustrated, to be less than 0.5 near ao — 0.4 
where there is little geometric frustration for WC forma- 
tion. Physically the reason that the overlap is large when 
the WC state is frustrated is that the Laughlin state at 
// ~ 1/7, while describing an incompressible liquid state, 
already has substantial short-range crystal order built 
into it; it is thus energetically still competitive. 
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FIG. 4. The square of the overlap between the Laughlin 
state and the exact ground state, in a system with aspect ratio 
ao = 3/8. It is plotted as a function of 6V5 (see text). The 
arrow on the right hand side indicates the position at which 
the near-degeneracy of the low-lying states starts to develop. 
The arrow on the left hand side indicates the position of SVs 
below which the absolute ground state is not in the K — 
sector. The overlap shown is between the Laughlin state and 
the lowest energy K = state. 

At = 1/7, the Laughlin state can be stabilized by 
slightly modifying the short range part of the Coulomb 
potential. In Fig. 4 we plot the square of the over- 
lap between the Laughlin state and the exact ground 
state of the system at ao = 3/8, as a funclian of SV^, 
the change of the m — 5 pseudo-potentialO from its 
Coulomb value. It is seen that the overlap is substan- 
tially improved to be above 0.8 by adding a small pseudo- 
potential SV^ = 0.025e^/e£ (which represents less than 
10% increase of the V5 pseudo-potential), even though 
the system geometry is most favorable for WC forma- 
tion. This is further indication that v = 1/7 is very close 
to the phase boundary separating the crystal and liquid 
phases. The smooth change of the overlap as a function 
of SVc, also suggests the transition is either continuous 
or very weakly first-order in the thermodynamic limit. 
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Adding a negative 5V5, on the other hand, further sup- 
presses the overlap; in particular, at 6V5 « — O.Olle^/e^ 
there is a level crossing for the ground state, and for 
SV5 < — O.Olle^/e^ the absolute ground state has a dif- 
ferent wavevector from the Laughlin state; thus the over- 
lap between the Laughlin state and the absolute ground 
state is strictly speaking zero here. However, the quasi- 
degeneraces are well developed and the ground state is 
not unique. We have tracked the overlap with the K — 
"ground state" . In real systems, the short-range part of 
the Coulomb interaction is softened due to the finite ex- 
tent of the electron wave function along the z direction 
(or finite layer thickness), thus the WC state is expected 
to be further stabilized; we have studied potentials of 
typical layer thickness and find it is indeed the case at 
I' = 1/7. 




FIG. 5. The real-space lowest Landau level projected 
(or guiding center) density-density correlation function, for 
N = 6 electrons at filling factor 1/7, in a rectangular system 
with aspect ratio 3/8. A 2D contour plot is also shown. 

As in our previous studies we now turn to the density 
response functions. In Fig. 5 we show the LLL projected 
(or guiding center) density-density correlation function, 
in real space. This is the Foueier transform of the guiding 
center static structure factortj. The central peak repre- 
sents the usual 6{r) singularity projected to the LLL, 
which is at the position of the electron where the cor- 
relation is measured from. The crystal structure can be 
clearly seen; the number of peaks (including the central 
one) equals the number of electrons iV = 6, indicating 
each unit cell includes one electron, which is what one 
expects for a Wigner Crystal. 

The tendency toward WC formation, as well as the 
structure of the resultant lattice, can also be seen in the 
density response function in momentum space, x(q)- Fig- 
6 is a plot of x(q) for the same system as in Fig. 5. Here 
we see strong response at an array of momenta, which 



form a lattice that is close to being triangular. This is 
essentially the reciprocal lattice of the crystal that forms 
in real space. The origin of the large response at these 
momenta is the near-degeneracy of the low-lying states 
that are connected by these momenta discussed earlier; 
small energy denominators lead to large response at these 
momenta. 
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FIG. 6. A 2D plot of the LLL projected density re- 
sponse function x(l)) with the size of the dot representing 
the strength. The system is the same as that if Fig. 5. The 4 
large dots represent a peak value of about 21000. The nearby 
smaller dots along the Qy axis represent a value of 4700 etc. 

In summary, we have shown that FQHE incompress- 
ible fluids become unstable around v = 1/7 in agreement 
with previous predictions. We find the Wigner crystal 
becomes the ground state as the system geometry is ad- 
justed to accommodate the crystal. It is more difficult to 
determine the order of the transition by the small sizes 
that we have studied. However, a strong level cross- 
ing first order transition seems to be ruled out by our 
study. We observe only a continuous evolution of the re- 
quired quasi-degenerate manifold which always includes 
the state a.t K ~ 0. This state for geometries that frus- 
trate the crystal shows large overlap with the Laughlin's 
state. This suggest the transition may be second order 
but a weakly fixst order transition predicted by a mean 
field approachllJ also remains a possibility. 

FDMH and EHR thank the Aspen Center for Physics 
where this work was started. EHR is grateful to NHMFL 
for their hospitality where it was brought to a conclu- 
sion. This work was supported by NSF grants DMR- 
9809483 (FDMH), DMR-0086191 (EHR), DMR-9971541 
(KY), and the A. P. Sloan Foundation (KY). 



4 



1 E. P. Wigner, Phys. Rev. 46, 1002 (1934). 
^ C. C. Grimes and G. Adams, Phys. Rev. Lett. 42 , 795 
(1979). 

^ For reviews, see, e.g., chapters by H. Fertig and M. 
Shayegan in Perspectives in Quantum Hall Effects, edited 
by S. Das Sarma and A. Pinczuk, Wiley, New York (1997). 

* H. Fukuyama, P. M. Platzman, and P. W. Anderson, Phys. 
Rev. B 19, 5211 (1979). 

® R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

® D. Yoshioka and H. Fukuyama, J. Phys. Soc. Jpn. 47, 394 

(1979); D. Yoshioka and P. A. Lee, Phys. Rev. B 27, 4986 

(1983). 

^ K. Maki and X. Zotos, Phys. Rev. B 28, 4349 (1983). 

* P. K. Lam and S. M. Girvin, Phys. Rev. B 30, 473 (1984). 
^ R. K. Kamilla and J. K. Jain, Phys. Rev. B 55, R13417 

(1997). 

° J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 

^ S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Phys. 

Rev. B 33, 2481 (1986). 
^ E. H. Rezayi, F. D. M. Haldane, and K. Yang, Phys. Rev. 

Lett. 83, 1219 (1999). 
^ F. D. M. Haldane, E. H. Rezayi, and K. Yang, Phys. Rev. 

Lett. 85, 5396 (2000). 

F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
^ F. D. M. Haldane, Phys. Rev. Lett. 55, 1095 (1985). 
F. D. M. Haldaue, in The Quantum Hall Effect, 2nd Ed., 
edited by R. E. Prange and S. M. Girvin (Springer, New 
York, 1990). 



5 



